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ABSTRACT 


A refined gravity field model, Goddard Earth Model 
GEM-L2, has been derived using the Lngeos orbital data 
yielding better baseline measurements for the analysis of 
tectonic plate motion. This field also contributes to an 
improved understanding of long wavelength features, such as 
the sea slope across broad ocean basins, through its 
significant improvement of the long wavelength geoid 
(through degree and order 4). The geoid for these terms has 
an accuracy estimated at 4^ 8 cm in GEM-L2. GKM-L2, as in 
all recent Goddard Earth Models, relies heavily on the 
precise near-Enrth satellite laser ranging data, in this 
case provided by NASA's Crustal Dynamics Program. Two and a 
half years of Lageos laser data acquired from over 20 well- 
distributed stations were combined with the existing data 
from the best satellite-derived model, GEM-9, to deveJ.op the 
new Lageos model. Testing shoves that the Lageos gravity 
field error at long wavelengths is less than half that for 
GEM-9. Independent tests using well determined longitude 
accelerations of 24-hour satellites have verified the 
Improved accuracy of the now model. A comparison of global 
laser "base" stations from independent data sets of 

alternating 15 day data segments over two years of Lageos 
show total inter-station positioning to ^ 1.8 cm when using 
this new field. The same comparison using the 1979 versus 
the 1980 Lageos data yields +_ 5.2 cm; this difference in 
agreement reflects the change in data distribution and other 
systematic errors along with the tectonic motion which 
has occurred between these chronologically distinct data 
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sets. Five day average polar motion values with a precision 
of 10 cm and change in length of day values accurate to 
better than .5 msec have been derived in the solution. The 
adjustment of these earth orientation parameters are 
necessary to achieve the accurate stations and geopotential 
results in GEM-L2. 
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1.0 INTRODUCTION 


The Lageos spacecraft is unique in several 
respects. In addition to its well-designed shape as a laser 
target, the satellite's high density and high altitude orbit 
virtually eliminate any error from the uncertainty arising 
in modeling non-gravltational forces. At an altitude of 
nearly one earth radius, Lageos experiences a strong 
gravitational force from the longest wavelength portion of 
the geopotential, with the effects falling off quickly at 
short wavelength due to the attenuation of the field at 
Lageos' altitude. 

Coupled with an extensive and highly accurate laser 
tracking campaign, Lageos data have provided a rich resource 
for Improving our knowledge of the long wavelength gravity 
field. Since the use of _Satellite J^aser _^anging (SLR) for 
measuring crustal dynamics requires the proper modeling of 
orbital dynamics, the gravity studies for Lageos have borne 
dual results: 1) an improved ability to estimate Intersite 

distances for tectonic motion studies, and 2) a greatly 
enhanced knowledge of the gravity field through degree and 
order four. 
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2.0 DATA 

The Lageos spacecraft was placed Into a high altitude 
(5900 kra) , nearly circular orbit on May 4 , 1976 , Typical 
mean Keplerian elements for Lageos are shown below; 

8eml~mnJor axis; 12267.7 kra 

eccentricity: .002 

inclination; 109?84 

To datej data from a global network of laser stations was 

acquired over a time span of six years on Lageos. Figure I 

presents a synopsis of the data available over the first 

four years (1976 through 1979) of Lageos tracking. As is 
apparent from this figure, there was a significant increase 
in data availability, Improved global distribution and many 
more participating NASA stations commencing in 1979. Of 
note, also, is the improvement of the SAO systems (systems 
with 7900 number designation) which began about 1979. 
Consequently, our gravity modeling activities utilized 

Lageos data beginning with February of 1979. 

The new model, GEM-L2, has incorporated 2 1/2 years 

of Lageos based data in combination v-;ith satellite tracking 
data taken on 30 other satellites as In GEM 9 (Lerch et al., 
1979). This solution, complete In spherical harmonics 
through degree and order 20, contained well over 600,000 
laser measurements, more than half of which were taken from 
Lageos. Lageos laser data are accurate down to about 8 cm 
and formed the predominant weight In the solution. However, 
because of the attenuation of the gravitational field at 


4 



ORIGINAL PAGE IS 
OF POOR QUALITY 


h- 

2: 

Hi 




0RIOINAL PA® IS 
ai: POOR QUAtnY 

high altitude the effective strongih of the Lagooa data 
falls off rapidly for the high degree terms and resulted in 
the GEM 9 data controlling the adjustment for the harmonics 
beyond degree 7. For Lageos approximately 20 laser stations 
with a global distribution contributed data in the GBM-L2 
solution although more than half of tliese stations were 
located in the United States. A description of the data 
( 1979 through 1981) is given by station for each of the 
orbital arcs (15 doys in length) in Tables la through Ic. 
The data on Lageos were weighted in the solution to provide 
for a balance both by station and by orbital arc as 
described in Table 2. 
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Tho potonUlnX coeEPtcients GKM-L2 

rofor to the tlePinltlon of the grivvltatlomii poCoivcloX as 

I 


who to 1 1 <|> , and X avo the distaiu'e to tho con ter of mass, 
ijeocentrle latitude and lonj’ltude. Is the fully 

normal, iKed associated hoRondre function of doRreo S, and 
order m, and r,^, is the earth's mean radius* Tho GKH-h?. 
model Is complete to degree and order 20, wltli Lageos being 
evaluated Cor its contrihutlon complete out through degree 
and order 16, and beyond for the aonals and certain other 
orders. Figure 2 shows the coefficients selected for 
Lagoos' goopotentlal contribution when combined with GEM- 
9. An estimate of the orbital perturbations on hagoos 
arising from the gravity field Is contained, also, In Figure 
2. The complete set of the GEM*'h2 coefficients are given in 
Table 3. 



Even though a large Lageos data set has been used in 
GEM-b2, this satellite only solution, like GEM-9 (herch et 
al., 1979), still requires the utilization of a modified 
least squares method to achieve stability in Its recovered 
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OBIGIMAL PA® ^ 

OF POOR QUAtlTY 

coofflclonts of tilsh degree. The raodified least squares 
method which is employed, is fully described in terch et 
al., 1977. An elaborate force model was utilized in the 
analysis of the Lagoos data. The solution is referred to a 
speed of light of 299792.5 m/sec because of the large data 
set in GEM-9 which used this pre-existing value. The 
Earth's M has been adjusted in the GEM-L2 solution 
(simultaneously with the gravity harmonics and station 
positions) yielding a value of 398600.607 Icra^/sec^.^ This 
value is in good agreement with other recent findings (e.g., 
Lerch et al., 1978). Solid earth tides were also modeled 
with (Love's numbers) h 2 * .60, A2“.075 and kg « .29. A 

solid earth phase angle of 2.018° was obtained in the 
solution . 

Also polar motion and Al-UTl variations were 
simultaneously estimated from 5 day segments of the Lageos 
data. A solar radiation pressure coefficient, Cj., was 
adjusted over each 15 day arc. The Jacchia (1971) 
atmospheric density model does not contain values at the 
altitude of Lageos and therefore an along track acceleration 
is also permitted to adjust in order to account for "drag- 
like" effects discussed by Rubincam (1980) and also Smith 
and Dunn (1980). 


^If the speed of light value of 299792.458 m/sec were to be 
used, the value of GM would be changed by -.167, equalling 
398600.440 km°/sec2. 
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4.0 EVALUATION OF THE RESULTS 


GEM-9 was the gravity model adopted initially In our 
LAGEOS analysis. GEM-9 (for the 1979 time frame, which 
coincided with the beginning of our effort) was considered 
the best available gravity model developed from a data set 
strictly composed of conventional satellite observations. 
GEM-9 did, however, use an augmented laser ranging 
observation set which included some of the newest NASA 
systems tracking GEOS-3 and Starlette. Therefore, GEM-9 was 
a larger and considerably improved field over preceding odd- 
numbered (which designate "satellite-only") GEM fields like 
that of GEM-7. Figure 3 presents the estimated uncertainty 
of the coefficients found in GEM-9. These error estimates 
were obtained from the scaled formal uncertainties of the 
solution obtained through an analysis against Independent 
data (see Lerch et al., 1975 and 1977). Based upon a 
calibration with surface gravity a scale factor of /To was 
derived to relate our formal noise only error estimates to 
the "true" uncertainty in the fields. Similar tests made 
with altimeter data confirmed that this scale factor was 
acceptable but somewhat conservative. The scaling factor 
was further confirmed from orbital data through the analysis 
of satellite resonant effects (Wagner and Lerch, 1978). 

Based upon this scale factor (/lO) the estimated 
accuracy of GEM-L2 is compared to GEM-9 in Table 4 for the 
longest wavelength terms in the model through degree and 
order 4. The actual coefficient differences between GEM-L2 
and GEM-9 are also shown for these same terms. The highly 
accurate laser tracking data, together with Lageos' unique 
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ORIGINAL FA^i !u 
OF POOR QUALITY 

aenoitivity to this portion of the field* Accounts for the 
Improvements indicated in Table 4. However, since GEM-L2 
contains the GEM-9 dato set, an independent vorlfleation of 
these accuracies was sought. An independent vorificntlon of 
these accuracies (mostly for terras where x-ra is even) has 
been performed at NGS/NOAA (Wagner, 1982). A comparison of 
24 hour satellite longitude accelerations observed through 
the study of mean elements was compared to those calculated 
from GEM-Ll (Lerch and Klosko, 1981) and GEH-L2 as shown in 
Table 5. The results for each of the eight 24 hour 

satellites are shown in Table 6 . 

The low degree and order portion of the gravity field 
cause long period changes in the longitude of deeply 
resonant (librating) 24 hour satellites. These 

perturbations commonly are seen as either long period 
libratlons of the longitude of the spacecraft from its 

nearly stationary position or as a secular drift in 

longitude. By evaluating several year histories of the 
spacecraft's mean Kepler coordinates and comparing these 
with the calculated trajectories from given geopotential 
models, one obtains a very good Independent measure of the 
accuracy of the fields at lowest degree and order 

where A-m is even thereby satisfying the resonance condition 
(principally terms CS2^2» ^ 83 ^ 3 ). Figure 4 shows 

the orbital evolution of the SYNC0M2 satellite (Wagner, 
1972) indicating that the magnitude of the effect of ^$ 2^2 
on the longitude of a stationary orbit can be several 
degrees over a few years. 
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TABLE 5 

UPDATED EVALUATION OF RECENT GRAVITY MODELS 


USING THE OBSERVED LONGITUDE ACCELERATIONS 
OF EIGHT 24-HOUR SATELLITES 


I 



SAO SE2 

GEM-9 

GRIM-3 

GEM-L1 

GEM-L2 


(1970) 

(1979) 

(1981) 

(1981) 

(1982) 

Weighted RMS Residual 
Accelerations 

15.7 

6.0 

21.6 

1,0 

.80 


Where the weighted residual aeeeleration Is (A_— A^<Ta )* 

u ^ 

* Aq Is observed longitude acceleration 

Aq is the calculated accelerations from gravity model 

is the estimated uncertainty ascribed to the observed longitude acceleration 

° (~1.x 10“8r/d2) 


CONTRIBUTION OF GEOPOTENTIAL BY DEGREE TO ACCELERATION 



CONTRIBUTION 

CONTRIBUTION 


FROM FULL 

FROM ESTIMATED 

i 

FIELD 

ERROR OF GEM-L2 

2 

3000X10“^ 

2.6 X 10“^ 

3 

400 

1.0 

4 

40 

0.2 

5 

4 

0.1 

6 

0.7 

0.0 


m 



TABLE 6 

gravity model evaluation using data from 24 HOUR SATELLITES 
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GEM“L2 overall satisfies these accelerations to their 
requisite accuracy. This is shown in Table 5 where the 

weighted RMS acceleration residual is equal to 1.0 for GEM- 
Ll and 0.80 for GEM-L2. The contribution of the 

geopotential by degree to the longitude accelerations are 
shown on the bottom of Table 5 for the full values of the 
harmonic coefficients and also for the errov estimates of 
Table 4. They show that while the accelerations are 
sensitive to the harmonic coefficients as high as degree 6, 
the GEM-L?. error estimates beyond degree 3 or 4 are not 
directly tested by these longitude accelerations. Hence, 
the accuracy of GEM-L2's higher degree terms cannot be 
tested by this method for they exceed the accuracy of the 
longitude accelerations themselves. However, since uo 
longitude acceleration residuals came close to approaching a 
2.6 sigma error for the degree 2 contribution (see Table 5), 
it was concluded that our stated accuracies for the 6x6 
portion of GEM-L2 are possibly conservative. 

The simultaneous estimation of polar motion and Al- 
UTl variations from the a priori 5-day mean values^ made an 
Important contribution to both the stability of the gravity 
solution as well as Improved station positioning. Separate 
gravity fields were made utilizing the 1979 and 1980 Lageos 
data, and were intercompared between themselves and with 
GEM-L2 . The new polar motion and Al-UTl values made 

an improvement of 30% in the RSS coefficient difference to 

^Bureau International De L'Heure (BIH) Circular D 5-day mean 
values using 90 day smoothing. Al-UTl =( Al-UTC)+( UTC-UTl ) 
where (Al-UTC) is from the U.S. Naval Observatory Bulletin 
and (UTC-UTl) is from BIH. 
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4x4 over fields holding the polar motion ond Al-UTl values 
fixed at 1)111 (e.g., the RSS coefficient difference between 
the 1979 vs. 1980 solutions decreased from 11 to 8 cm with 
the adjustment of polar motion ~ see Table 7). The 5-day 
mean polar motion values themselves were well determined in 
GEM-L2 and have a precision of about 10 cm. Figures 5 and 6 
show examples of the corrections made to BIH (Circular 0 
values) polar motion in GEM-L2 for the X and Y coordinates 
during 1980. Adjustments of up to 96 cm are seen from the a 
priori values. The full set of GEM-L2 polar motion 
and ALOD (changes in length of day) computed from the Al- 
UTl values are found in Tables 8 and 9. 

In the procedure with laser data, even after one 
station's longitude is hold fixed, a singularity still 
exists for each orbital arc between the positioning of the 
satellite's longitude (node) and tlte earth's rotation 
(Al-UTl). The three 5-day mean values of Al-UTl in each 
15-day arc are good to within an undetorminod constant and 
henc^ one value is held fixed at the a priori value. The 
three mean values are differenced to provide two accurate 
values of 5-day mean ALOD and are listed for each arc in 
Tables 8 and 9. In order to benefit from the adjusted Al- 
UTl values in the computation of orbits, the specific arcs 
used in GEM-L2 or arcs contained v^ithin these 15-day 
intervals must be employed due to the discontinuity. Using 
this procedure the values recovered in GEM-L2 were used in 
our orbital analysis and provided the means for correcting 
the a priori values of Al-UTl in our test arcs. The use of 
this corrected information greatly improved our orbital 
results as subsequently shown. 
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FIGURE 6 
POLAR MOTION: 
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Satellite positioning accuracy directly reflects on 
the capabilities of using SLR to monitor tectonic motions 
and deformations. GEM-L2 , by improving the force modeling 
on the Lageos orbit, yields a more accurate epheraeris than 
previously available, which in turn Improves the accuracy of 
the estimated station coordinates and their inter~station 
distances. The Improvements in orbital accuracy are 
discussed at length in the next section. Typically, the 
post-fit RMS of the range residuals (normal points) on 
Lageos when using GEM-L2 in a 15 day arc is on the order of 
10 to 15 cm. Consequently polar motion and Al-UTl errors, 
which produce range residuals of this magnitude or in excess 
of it, map into the positioning of the stations with respect 
to the orbit and must be accounted for if stable station 
position solutions are to result. This is clearly shown in 
the results presented in Table 10. 

Table 10 shows a comparison of laser interstation 
distances for 8 "base" stations using two Independent solu- 
tions. The values in Table 10 are the ave7.-age error in the 
magnitude of the baseline differences. The Lageos data set 
(from 1979 through December of 1980) has been divided into 
two parts where each single month's data is split into two 
15 day segments and each alternating segment contributes to 
one of these two independent solutions. Therefore, while 
the two solutions contain independent data, they both span 
the same time Interval, have nearly equal representation for 
each station, and should average out (for the purposes of 
this comparison) any resulting plate motion over these 
baselines. The impact of the GEM-9 and GEM-L2 force models 
on station positioning are intercompa red as well as the 
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TABLE 8 

PriL,RR MOT I OH FROM GFM-LS 


POCH TIME 

EHD TIME 

y. OF POLE 

y OF POLE 

(VVMMDn.l 

(VYMMriD) 

mpc SEc:i 

(RRC 6EC:i 

?9eisen 

790206 

fi. 04236 

0.06581 

79QB36 

790211 

0, 03006 

0.07124 

790ai 1 

790216 

0.00962 

U. 06*1*23 

790216 

790221 

-.01840 

0.07566 

790221 

790226 

0.00731 

0. 06682 

790226 

790303 

-.00234 

0 . 06454 

790303 

79030.8 

• 0ill4izl f’ 

0. 08268 

790303 

790313 

-.04352 

0. 10821 

790313 

790318 

. 03390 

0. 10377 

7903 1 8 

730333 

-*.07209 

0. 097:30 

7903S3 

790328 

-. 10478 

0. 12547 

790328 

790402 

10689 

0. 1375? 

790402 

790407 

-.11646 

0 1 1 420 1 

790407 

790412 

~. 12197 

0. 16053 

790412 

T’90417 

-. 12915 

0. 16711 

790417 

790422 

-. 13187 

0. 17149 

7904cli::! 

790427 

-. 13487 

0.20258 

790427 

790502 

-. 1:3713 

0.20528 

790502 

790507 

-. 14700 

0. 23164 

790507 

790512 

• 1 4cI3cI 

0. £3456 

73bl5i ci 

790517 

-. 14868 

0 . 23378 

790517 

7’90522 

-. 15254 

0.25760 

790522 

790527 

-, 15193 

8 . 26936 

790527 

79060 1 

14911 

0 . 27508 

790603 

790608 

-. 14744 

0. 29402 

790608 

7906 1 3 

— . j '3 2 6 -5 

0.28258 

790613 

790618 

-. 13930 

0. 31097 

790618 

79062:3 

0. l2Wt'4* 

0.24472* 

790623 

790628 

-. 12911* 

0. 34322* 

790628 

790703 

-.14598* 

0.41153* 

790703 

790708 

-. 09748 

0.:34779 

790703 

790713 

-.09682 

0.:36196 

790713 

790718 

-. 10153 

0.:3?976 

790718 

790723 

- . 08346 

0.:38751 

790723 

790728 

— . 008 1 5 

0. 38973 

790728 

790802 

06076 

0. 39404 

790802 

7'3L10U7' 

-.04733 

8 .41 488 

790807 

7908 1 2 

-.02678 

0. 42905 

790812 

790817 

-.01391 

0.4.2019 

790817 

790822 

-.00727 

0. 41940 

790822 

790827 

0 . 0 1 385 

8 . 42778 

790827 

790901 

0.01165* 

0. 42513* 

790901 


0.07187 

0.41025 

790906 

790911 

0.05219 

0. 43264 

790911 

790916 

0.05721 

0 . 4..T'2c'8 

79y 9 1 6 

790921 

0 . 07397 

8 . 42:3:35 

790921 

790926 

0.07021 

0.43368 

790926 

791001 

0.09166 

0. 42970 

791001 

79 1 006 

0 . 09349 

0. 40731 

79 1 006 

791011 

0. 10094 

0.41176 

791011 

791016 

0. 11009 

0.:59841 • 


Insufficient data; poorly determined values 
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TABLE a (CONTINUED) 

POLhP tlCiTJOH FROM GEf1-lc‘ 


ORIGINAI- Pf|Sf ® 
OF POOR QUALITY 


EPOCH TIME 

ENP time 

OF POLE 

OF POLL 

iVVHMnri' 

'VVMMDPj 

iflF'C CEC) 

1 Hpr 8EC 1 

r-jiciic 

7-'* 1021 

0.11346 

0. 388^31 


78102t. 

0.11831 

0 . 38255 

roiosc 

781051 

0. 13151 

0. 3031 I 

roio'ii 

781105 

0. 13802 

»*’. 37152 

TO 1105 

781110 

0, 14515 

0. 36058 

roll 10 

781115 

0. 14818 

0. OOOT-^ 

701115 

781120 

0. 15053 

0.34341 

701 lao 

781125 

0. 13841 

0. 33747 

701 las 

781130 

0. 14077 

0.32170 

701130 

781205 

0. 14870 

0. 31184 

701205 

?-*n210 

0, 14685 

0. 30421 

7012)0 

781215 

0. 14422 

0. 28*476 

701215 

781220 

0. 1449i2 

0. 28006 

"’01 220 

781225 

0. 14258 

0.27804 

701225 

783230 

0. 15174 

U » 

800102 

1500 1 07 

0. 15036 

0.£'5186 

800107 

2001 12 

U « i 444u 

0.24 ■’06 

8001 12 

800117 

0. 13140 

0.24356 

800117 

yiJU 1 c!cl 

0. 12183 


800122 

1 £7 


0 . u26t5 

800127 

800201 

0. 11801 

0. 220-6 

80020 1 

S00c:06 

0. 1 1028 

0. 20601 

s“JOO£0b 

80023 1 

0.08687 

0 . 20555 

0002 1 1 

800216 

0.08302 

0 . 20381:! 

800216 

800221 

0* 07b6? 

C6 18528 

800221 

SOO£c*b 

0 . 07467 

0. 18268 

8U02c't' 

SCifi 

0.07474 

0. 18840 

800302 


0.07374 

0. 10160 

800307 

oori:ii£ 

0.06761 

0 . 1 8773 

800312 

0OC1717 

0.05152 

0 . 1 836? 

800317 

sooy52 

0.04221 

0. 18745 

800322 

J:Ci03£r 

0.02451 

0. 18884 

000327 

000401 

0.02357 

0. 18848 

80040 1 

800406 

0. 02056 

0. 18245 

800406 

800411 

0.00381 

0. 20385 

800411 

8004 1 6 

- . 00724 

0 . 20856* 

800416 

800421 

-.08151 

0 . 20750 

800421 

800426 

-.02221 

0. 2176*8 

800426 

800501 

no*5»oo 

« UC 1 s«'W 

61 . 22346 

800501 

800506 

— . 028fc'5 

0 . 22351 

800506 

80051 1 

-.03824 

0.23284 

80051 1 

£!Ci05 1 6 

-.04372 

0.23680 

800516 

W0L^52 1 

— . Oc'84 1 

0.23761 

800521 

800526 

— • U4502 

0 . 2448 1 

800526 

800531 

- . 04830 

0.25286 

800531 

800605 

-.04184 

0.25608 

800605 

8006 1 0 

^ r*4Cro*T’‘5’ 
”” « Uv»C. 1 1 

0.26201 

800610 

800615 

- . 04652 

0. 26776 

800615 

800620 

-.05282 

0.27515 

800620 

800625 

- . 046 1 2 

0.276*12 

800625 

800630 

- . 0383 1 

0 . 28083 
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POLFIP HOT I OH FROM GEM-LS 


EPOCH TIME 

Ef(n TIME 

M OF POLE 

V OF POLE 

(VVMMUrO 

iyvmmdd:i 

I'FIRC SEC:i 

tpRC SECi 


800701 


0551 1 

0.29827 

800706 

800711 

• U4 1 hui 

0 i £8676 

©00711 

800716 

03260 

0.30141 

800716 

800721 

- . 03657 

0.30816 

300721 

800726 

-.03092 

0.31439 

©00726 

800731 

- • 033SS 

0.32159 

800731 

800805 

- . 03553 

0.31654 

800805 

8003 10 

02307 

0. 32481 

800810 

800815 

- . 03024 

0.32520 

8008 1 5 

800820 

-.02572 

0.32693 

800820 

800825 

-.01953 

9 • 334 S3 

800825 

800330 

■ yC'JC'O 

0. 33547 

800830 

800904 

-.02577 

0.34264 

800804 

800809 

-.01982 

0. 34304 

800808 

800914 

-.02404 

0.34017 

800814 

800919 

-.01789 

0.34540 

800818 

Si00924 

-.01904 

0.35156 

800824 

800929 

-.02451 

0 . 35282 

800828 

80 1 004 

90803 

0.35642 

801004 

80) 009 

-.01127 

0.35375 

801008 

801014 

-.01399 

0 . 35480 

801014 

801019 

- . 0068 1 

0 . 37008 

80 1018 

801024 

-.00160 

0.36526 

801024 

80 1 029 

-.00945 

0.37419 

801 028 

801103 

0.00243 

0.37585 

801 103 

801108 

0 . 00587 

0 . 36965 

801 108 

801113 

0. 01022 

0.38007 

801 113 

801 US 

0.01466 

0.37843 

801 1 IS 

801123 

0.02706 

0. 37407 

801 123 

yo 1128 

0. 03433 

0 . 38633 

801 128 

301203 

0.04107 

0.37814 

801203 

801208 

0 . 044 6* 3 

0.37785 

801 208 

801213 

0 . 05022 

0 . 38555 

801213 

801218 

0.05987 

0.37261 

801218 

801223 

0.06229 

0.37020 

801223 

1 228 

0.07107 

0. 37347 

801 228 

810102 

0. 07320 

0 . 36536 

810102 

810107 

0 . 07508 

0.36159 

810107 

810112 

0 . 07839 

0.37329 

810112 

810117 

0.08586 

0 . 35664 

810117 

810122 

0 . 07958 

0.35986 

810122 

810127 

0 . 08426* 

IJ » JOD i* 5 

810127 

810201 

0. 08829 

0 . 34396 

810201 

810206 

0.09155 

0.33750 

810206 

810211 

0.09535 

0.33137 

81021 1 

810216 

0.09981 

0 . 33657 

810216 

y 1 0221 

0 . 09405 

0.33510 

810221 

810226 

0.09683 

0.31882 

810226 

81 0303 

0. 10090 

0. 31807 

810303 

81 0308 

0.08837 

0.32041 

810308 

810313 

a 

0. 09328 

0 . 30899 
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POLAR HOT I on FPOH GEH~La 


EPOCH TIME 

EHD TIME 

M OF POLE 

V OF F'CiLE 

iVVMMriD'i 

iVVMMUD^ 

(ARC SEC’i 

t APC SEIM 


SI 03 13 

810318 

O.O8?80 

0.31154 

SI 03 18 

8 10323 

0. 10133 

0.30774 

8103S3 

8 1 0328 

0. 11215 

0. 28541 

810320 

8104132 

0 . 08384 

0 . 30843 

SI 0402 

81040? 

0.08140 

0.28810 

81040? 

810412 

0. 10261 

0. 28538 

SI 04 12 

01041? 

0.0808? 

0.28124 

31041? 

810422 

0. 08544 

0.27851 

810422 

81042? 

0. 1064? 

U . clt’4'"i3 

S1042? 

810502 

0. 0854? 

0. 27783 

810502 

81050? 

0. 10013 

0 a C*6‘ 1 64 

81 050? 

810512 

0.08482 

0.25310 

010512 

81051? 

0.08274 

0 . 26566 

81051? 

810522 

0. 10048 

0.25008 

810522 

81052? 

0. 0837? 

0.24716 

81052? 

810601 

0.1124? 

0.24830 

810601 

ij 1 06'0t' 

0. 10014 

0.24548 

S 1 0606 

810611 

0 . 0870? 

0. 23571 
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TABLE 0 

CHl'INGE IN LENGTH OF Dfl’i' FROM GEM-LE 


MIDPOINT DATE 

DELTA LOn 

(YVMMDIO 

i:5EC:i 


?90£G6 

. O0S80 

7O0211 

. 00S76 

?00££i 

. 00306* 

?'?0££6 

. 00360 * 

?'?03O8 

. 0O£83 

700313 

. 00£8 1 

70O3£3 

. 00356 

rvoco 

. 00370 

70040? 

. O0£7? 

7004 1£: 

.00316 

7004S2 

. 00374 

7004S7 

. 00330 

700Si07 

. OO£00 

7005 1£ 

. 0O£08 

?'005££ 

. 00300 

7005£7 

. 00£50 

700608 

.00841 

7008 1 3 

. 00££iJ:* 

7006£3 

. 00108# 

7008££ 

. O0£0£# 

700708 

. O0££ 1 

700713 

. 00£46 

700? £3 

.00151 

70G7£8 

. 00182 

70080? 

. O0£4S 

70uy 1 fci 

. 00266 

7008££ 

. 00100# 

?008£7 

. O0£O0 

700006 

. O0232 

7000 1 1 

. O0£60 

?0O0£1 

. 00247 

?000£6 

. 0O££0 

701006 

. 00321 

701011 

. 00264 
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IiPOIHT IiRTF 

ivvNMiiri:i 


riELTft Lcm 

(SEc;i 




791105 

79 11 £5 

?91£05 
79 IS 10 

791S£0 

791SS5 

900107 
9001 IS 

900 1££ 
9001 £7 

900£09 

900S11 

900££1 

900Sc'9 

800307 
8003 IS 

8003££ 

O003S7 

800406 
8004 1 1 

8004S1 
8004 £6 

800506 
8005 1 1 

8005S1 
80 05 £6 

8006015 

800610 


W 006 c:w 

8006£5 


’ . 00S80 
. 00£58 

. 00£69 
.00S14 

.00S51 
. 00S46 

. 00S35 
. 00S49 

• l3U£43 
. 00£4£ 

. 00£47 

• 010S49 

. 0OS89 
,00£15 

. O0£££ 

.00199 

. 00£7£ 

. 01 01 S cl £ 

.00SS1 
> L10SS 1 

. 00£60 
. 0084 1 

. 00S4S 
. 0OS70 

. 00££9 
. 0i0i£55 

. 00S38 

. 00S86 

. 010S15 
. 010S14 

. O0£14 
. 00£19 

• 010 176 
. 00160 
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1 


mmc lu 

LEMGTH OF P0‘i' PODH 

8EH-L0 




MIIIPDIHT IiFiTE 


X'Fl.TO LOU 




ivvMMrmi 


1 OEi: f 




S’Dinrcifc. 


.8810l- 


! 

i 


pr.iurii 


,U81J? 


, 




. 8818 !! 


1 




. 88130 




!"iLiCiCiD'? 


, 88181 




iBCiD©lG 


. i:i8 1 0£ 


1 


?:C;iCi5{£L‘D 


. 88188 


fi 

1 


P00825 


. 8i!t 1 08 




800884 


. 811100 




mm9 


. 88£1? 




080818 


,88£10 


\ 


S88884 


. CinOOD 


''4 


081084 


.88018 


■J 

1 


081888 


. 88e'r>0 


1 


881818 


. !i8£i;.8 


I 

\ 


l•:!818£4 


. 88 ;i££ 


} 

j 


881 18H 


> 88£t.'£ 


! 


I'lO 1 1 88 


. 88808 


1 


081 1 18 


. 880T4 




881100 


. 88081 




881 £00 


. 88£40 




081080 


. 8801 1 




081 £10 


. 88000 




i.nO 1 w£'J 


. 88£1 4 




01818E 


. 88££4 
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TABLE 9 (CONTINUED) 


;:HFlHGE IN LENGTH OF DRY FROM GEM-L2 


MinPOIMT DRTE 
i:vYMMrin:i 


DELTA LOD 
(SEC.1 


810318 

010333 

8 1 04 0S 
81040? 

81041? 
8 1 04cIir.I 

810503 
8 1 050? 

81051 ? 

810533 

810801 
8 1 0606 


. 008S5 
. 0035? 

. 0030 1 
. 0030 1 

. 00310 
. 00S6S 

. 00334 
. 00254 

. 00245 
. 00228 

. 00295 
.00218 
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a priori modeling of polar motion and Al-UTl with that 
obtained from LAGEOS for their contribution to station 
positioning. An excellent 1.8 cm agreement is obtained for 
the baseline differences of these two solutions when using 
the gravity and Lageos earth orientation parameters of the 
GEM-L2 solution. The 28 baseline-comparison from the GEM-L2 
model is shown as a histogram in Figure 7.^ All 28 baseline 
differences from these alternating 15-day arcs are shown in 
Figure 8. GEM-9 with a priori polar motion and Al-UTl 
values, which reflects the state of the art prior to the 
Lageos analysis as performed by ourselves and many others, 
yields about a four-fold worse agreement. These results 
clearly show that the new LAGEOS polar motion and Al-UTl and 
the gravity of GEM-L2 are needed to achieve the 1.8 cm 
baseline comparison. 

A comparison was also made dividing the data 
chronologically into a 1979 versus a 1980 set of 
solutions. In an ideal environment, the sole difference 
between this comparison and that described in Table 10 would 
be the plate motion which has resulted over this year. 
However, In this later case, data distribution and station 
participation also varied and must be considered. However, 
the results are still quite satisfying and are shown in 
Table 11. A 6 cm average has been obtained for the total 
position differences over these two years when using GEM- 
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^The laser system at Westford, Mass, (station number 7091) • 1 

seems to be somewhat inconsistent in longitude as compared J 

with the other sites. The longitude coordinate is most 
effected by any range bias. See Chr is t odoul id is and Smith, 

1981 . 1 
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36 


I 




ORIGINftU PftGE ra 

OF POOR QUAUTY 


L2 . With GIiM~L2 the b«8oline ecabllity over the Uivltod 
States is 4tl cm between 1*>70 and 1980 solutions, afiain, a 
good fcsult. 

An ottor analysis has been performed to give an 
estimate of the expected baseline error due to the estimated 
uncertainties in the GEK-h2 model for tlioso 28 "base 
station" baselines. In this analysis, the actual normal 
matrix obtained from 2 years of LAGEOS tracking was 
employed. Table 12 presents these resvilts indicating that 
the gravity field error in GEM-1,2 Is the cause of 1.7 cm 
average error on these baselines. However, 2 years of data 
is averaged (in a geometrical sense about these stations) to 
form these results. Also the gravity field error should 
become somewhat worse for more limited data distributions 
found in annual or monthly solutions. 

The tracking station coordinates obtained from GEM-L2 
are presented in Table 13. Those coordinates refer to the 
optical axes of the laser instruments. Tltese station 

coordinates, unlike the ones discussed above, represent 
values derived from the. entire 2 1/2 years of hageos 
tracking. As such, they most likely average out any 
interceding plate motion over 1979 to 1981. They, however, 
represent a starting set of consistent coordinates to be 
utilized witli GEH™li2 . 
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FIGURE 7 

HISTOGRAM OF BASELINE AGREEMENT FOR 8 BASE STATIONS 
(AVG - 1 8 cm) WHEN USING GEM-L7 GRAVITY AND POLAR MOTION 
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FIGURES 

A MEASURE OF SLR BASELINE PRECISION 
FROM ALTERNATING ARCS OF 
1979 THROUGH 1980 LAGEOS DATA 


i! 




STATION NO. 



7086 

7090 

7091 

7114 

7115 

7907 

CO 

r* 

7063 

2.2 

-0.5 

3.0 

2.6 

-0.5 

-2,4 

-0.4 

7086 


3,6 

4.9 

2.0 

-1.0 

-1.3 

3.3 

7090 



-0.7 

1,1 

1.0 

-0.0 

-3.7 

7091 




4.5 

1.6 

-O.S 

-0.1 

7114 





1.4 

-1.0 

0.3 

7115 






-3,5 

0,8 

7907 







2.3 


ACTUAL BASELINE 
DIFFERENCES (CM) 


i 


BASELINE AGREEMENT FOR 8 BASE STATIONS 
(AVG = 1,8 cm) WHEN USING GEM-L2 GRAVITY AND 
LAGEOS POLAR MOTION 
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TABLE 11 

GEM-9 VERSUS GEM-L2 STATION POSITIONING COMPARISON 
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AND A1-UT1 FROM RESPECTIVE SOURCES. BIN CIRCULAR D 90-DAY SMOOTH VALUES 


TABLE 12 

BASELINE UNCERTAINTIES COMPUTED FROM 
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TABLE 13. GEM-L2 STATION COORDINATES 
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5.0 ORBITAL ACCURACIES ON LAGEOS 

Dynamic satellite geodesy implies that improvement in 
the modeling of the forces acting upon the satellite should 
yield less error in the estimation of station positions and 
parameters describing the location of the earth with respect 
to the satellite's orbital plane. With laser systems 
routinely producing range observations at the 10 cm level of 
accuracy, a great deal of improvement was possible in our 
modeling of the geopotential to fully exploit this data to 
improve the capability of measuring plate tectonic 

activity. This was the principal direction of our work. 
Initial error analysis mapping the uncertainties found in 
each of the coefficients of GEM-9 revealed that LAGEOS 
orbital errors of ± Im or more could be expected. This 
analysis was made using the linear perturbation theory 
developed by Kaula (1966, pp 40 and 49). Each coefficient's 
uncertainty was used to estimate the magnitude of the 
resulting total perturbation error on LAGEOS for each 

harmonic constituent. Since most of the simple spherical 
harmonics give rise to a number of perturbations at 

different frequencies on LAGEOS' orbit, these were combined 
as their RSS error for each term. Figure 9 shows the 
results of this analysis for the GEM-9 uncertainties. In 
Figure 9, the low degree and order field comes through as 

the strongest source of error witli the exception of the 
resonance terms of order 6 and 7 (corresponding to LAGEOS's 
6.4 orbital revolutions per day). The total RSS of all of 
the errors shown in Figure 9 is approximately ± Im. This 
analysis was repeated at tlie conclusion of the development 
of GEM-L2 and is shown for the estimates of the GEM-L2 


44 


FIGURE 9 

LAGEOS ORBIT ERRORS (CENTIMETERS) 
DUE TO ERRORS IN GEM 9 HARMONICS 


ORIGINAU PAGE 13 
OF POOR QUALITY 



-Uti ^ ^ 


OTOGINAI. PAG 
OF POOR QUAUTY 


coefficient errors in Figure 10. Based upon this analysis, 
LAGEOS orbital accuracies lie In the vicinity of 35 cm when 
using GEH-L2. 

A direct measure of the orbital Improvement gained 
through the utilization of GEH-L2 has been obtained through 
a series of tests using actual LAGEOS data. Some of these 
tests are very severe they are designed to give a 
magnified insight into orbit errors arising from the gravity 
field. Therefore, the accuracies which are seen may be 
considered as worse case phenomena, for they are tests where 
large data gaps have been imposed on the Lageos observation 
set. In the tests shown in Table 14, a fifteen day span of 
LAGEOS data is selected. The middle five days of data is 
then deleted so that tlie data which is utilized for the 
original orbit computations consist solely of the data 
available from the first and last five day data spans. This 
fifteen day test orbit over this middle five day interval is 
therefore completely independent of the data taken during 
this interval. This independent five day span of data is 
then used to calculate a second orbit. This five day orbit 
is differenced over its interval with that of the original 
fifteen day arc in three components — radially, along track 
and across track. Table 14 shows a sample (representative) 
of the results we have obtained. Over the two year period 
tested, GEM-L2 outperforms GEM-9 in all cases and gives 
results averaging 26 cm in total orbit position error which 
closely reflects the magnitude of error expected from our 
error analysis using the linear perturbation theory. These 
improved orbital tests of necessity Incorporated the GEM-L2 
polar motion and Al-UTl values since the a priori values 
were corrected for these 15 day arcs. It is this 


FIGURE 10 

LAGEOS ORBIT ERRORS (CENTIPIETERS) 
DUE TO ERRORS IN GEM-L2 HARMONICS 
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Improvomoiit In the. modeling of the grnvitationnl forcee on 
liAGBOS which onhancea the ntilinntion of LAGBOS Inner data 
for Ulve estimation of the tectonic plate motion and dynamics 
of the earth/ s rotation and polar wandering • 

A second series of runs was made showing (Table 15) 
the RHS of fit to "normal" point laser ohservatlons created 
from the original range observation set. The use of GEM-L2 
and the LAGGOS polar motion/Al-UTl shows a factor of two 
improvement in our ability to fit these nearly "noiseless" 
observations. A typical RMS of fit to the normal points in 
15 day ares is about 10 cm when using GRM-h2 with LAGGOS 
polar motion and Al-UTl. 
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TABLE 14 

LAGEOS ORBIT TEST: 

!■ I 

16 DAY ARC WITH MIDDLE 5 DAY'S DATA DELETED |i J 


ORBIT FROM MIDDLE 5 DAY 


RMS SATELLITE POSITION DIFFERENCES 


EPOCH 

GRAVITY 

RADIAL 

ACROSS 

TRACK 

ALONG 

TRACK 

TOTAL 

DATE 

MODEL 

(cm) 

(cm) 

(cm) 

(cm) 

3/80 

GEM-9 

7.9 

81.7 

59.7 

101.5 


GEM-L2 

3.4 

17.9 

27.1 

32.7 

2/80 

GEM-9 

20.4 

78.0 

270.0 

281.8 


GEM-L2 

2.2 

7.6 

13.1 

'•5.3 

12/79 

GEM-9 

17.6 

15.2 

73.0 

76.6 


GEM-L2 

10.7 

4.6 

23.0 

25.8 

11/79 

GEM-9 

5.1 

31.8 

82.6 

88.7 


GEM-L2 

10.7 

32.9 

41.8 

54,3 

10/79 

GEM-9 

4.6 

61.1 

205.0 

211.3 


GEM-L2 

5.4 

31.9 

41.9 

52.9 


Estimated error in each orbit is 26 cm. As an example, since the error of 32.7 cm for 
the first case reflects errors from two orbits, the estimated error for each of the orbits 
(assuming they are equal) would bo 32.7 -r-v/2. 
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TABLE 15 

15 DAY RMS OF FIT TO NORMAL POINT OBSERVATIONS ON LAGEOS 


POLAR 


EPOCH 

GRAVITY 

MODEL 

MOTION 

A1-UT1 

NO. OBS. 

RMS OF 
FIT (CM) 

791215 

GEM-9 

BIH* 

602 

17.3 


GEM-L2 

LAGEOS 

602 

8.0 

790901 

GEM-9 

BIH 

203 

12.2 


GEM-L2 

LAGEOS 

203 

6.6 

800416 

GEM-9 

BIH 

663 

33.6 


GEM-L2 

LAGEOS 

663 

16.7 

800716 

GEM-9 

BIH 

771 

39.5 


GEM-L2 

LAGEOS 

771 

11.3 

800830 

GEM-9 

BIH 

1376 

36.7 


GEM-L2 

LAGEOS 

1376 

13.8 


*BIH Circular D 90— day smooth values was used as a priori. 
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6.0 SUMMARY 

A new gravity model GEM-L2 has been derived as part 
of the analysis being undertaken on Lageos for the NASA 
Crustal Dynamics Program. Lageos' unique contribution 
toward the resolution of the long wavelength (for terms 
through degree and order 4) geopotential has resulted in a 
two-fold Improvement in our knowledge of these terms. This 
Improvement is confirmed through an analysis of 24-hour 
synchronous satellite orbits. Over 400,000 Lageos ranges 
taken by 20 laser tracking stations over a 2 1/2 year period 
have been utilized in the new field. GEM-L2 contains a 
consistent set of parameters Including a geopotential field 
(20 X 20), tracking station coordinates (GEM 9 and 
Lageos), 5-day polar motion and ALOD values, and GM = 
398600.607 km^/sec^ based upon c = 299792.5 m/sec. 

GEM-L2 provides for better baseline recovery, orbital 
modeling, and station positioning from Lageos tracking. The 
adjustment of earth orientation parameters was essential for 
improving the geopotential and producing better interstation 
distances. The accuracy of the station positioning results 
(4^ 6 cm in total position globally and +.1*8 cm in 
baselines) of this report confirms the utility of using SLR 
for the measurement of tectonic plate motion and inter-plate 
deformations. These results are summarized in Table 16. 
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SUMMARY OF GEM-L2 
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